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ABSTRACT 


To  a  good  first  approximation  acoustic  propagation  in  an  underwater  sound 
channel  is  dominated  by  a  finite  number  of  trapped  modes.  However,  exact  solutions 
are  known  for  only  a  few  special  cases,  making  it  necessary  in  general  to  use 
numerical  methods  to  solve  the  normal  mode  equation.  But  often  one  is  interested  only 
in  the  gross  features,  such  as  the  number  of  modes  or  cutoff  frequencies,  and  one 
does  not  need  the  detail  provided  by  a  complete  normal  mode  calculation.  Even  if  a 
normal  mode  calculation  is  desired,  the  computation  time  can  be  reduced  considerably  if 
the  mode  wavenumbers  can  be  estimated  in  advance.  In  such  a  case,  the  WKB  method 
can  be  used  to  obtain  formulae  which,  although  they  are  approximate,  are  given  in 
closed  form.  In  this  paper  formulae  based  on  exact  and  WKB  solutions  are  presented 
for  the  number  of  modes  trapped  in  some  simple  sound  channels  and  for  the  wave 
numbers  and  cutoff  frequencies  associated  with  these  modes.  The  number  of  trapped 
modes  is  shown  to  depend  on  the  gross  features  of  the  sound  channel,  while  the 
distribution  of  modal  wave  numbers  depends  to  a  greater  degree  on  the  details  of  the 
sound  speed  profile  shape. 
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RESUME 


Oans  une  bonne  premiere  approximation,  la  propagation  acoustique  dans  un  canal 
de  son  sous-marin  est  dominee  par  un  nombre  fini  de  modes  pieges.  Cepandant,  des 
solutions  exactes  ne  sont  connues  que  pour  quelques  cas  speciaux,  ce  qui  oblige  en 
general  a  utiiiser  des  methodes  numeriques  pour  resoudre  I'equation  du  mode  normal. 
Souvent  pourtant,  le  chercheur  ne  s'interesse  qu'aux  caracteristiques  brutes  comme  le 
nombre  de  modes  ou  les  frequences  de  coupure  et  il  n'a  pas  besoin  de  la  quantite  de 
details  fournie  par  un  calcul  complet  du  mode  normal.  M@me,  lorsqu'un  calcu  du  mode 
normal  est  vouiu,  le  temps  de  calcul  peut  etre  considerablement  reduit  si  les  nombres 
d'onde  du  mode  peuvent  Stre  estimes  auparavant.  Dans  un  tel  cas,  la  methods  WKB 
peut  servlr  a  obtenir  des  formules  qui,  bien  qu'elles  soient  approximatives,  se 
presentent  sous  une  forme  fermee.  Dans  cette  communication,  I'auteur  presente  des 
formules  basees  sur  des  solutions  exactes  et  sur  des  approximations  WKB  pour  le 
nombre  de  modes  pieges  dans  des  canaux  de  son  simple  et  pour  les  nombres  d'onde  et 
les  frequences  de  coupure  lies  a  ces  modes.  II  demontre  que  le  nombre  de  modes 
pieges  depend  des  caracteristiques  brutes  du  canal  de  son,  tandis  que  la  distribution 
des  nombres  d'onde  modaux  depend  dans  une  plus  grande  mesure  des  details  de  la 
forme  du  profil  de  la  vitesse  du  son. 
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1  Introduction 


Acoustic  propagation  in  an  undorwatar  sound  channel  is  dominated  by  a  finite  number  of 
trapped  modes  whose  wavenumbers  depend  on  the  sound  speed  profile  in  the  channel.  Exact 
solutions  are  known  for  only  a  few  special  profiles,  making  it  necessary  in  genaral  to  use  numerical 
methods  to  solve  the  normal  mode  equation.  But  often  one  is  interested  only  in  the  gross  features, 
such  as  the  number  of  modes  or  the  cutoff  frequencies,  and  one  does  not  need  the  detail  provided 
by  a  complete  normal  mode  calculation.  Moreover,  even  if  a  normal  mode  calculation  is  desired,  the 
computation  time  can  be  reduced  considerably  if  the  mode  numbers  can  be  estimated  in  advance. 
In  such  cases  the  WKB  (after  Wentzel-Kramers-Brtllouin  and  others)  method  can  be  used  to  obtain 
formulae  which,  although  they  are  approximate1,  are  given  in  closed  form. 

In  this  paper,  formulae  based  on  WKB  solutions  are  presented  for  the  number  of  modes 
trapped  in  some  simple  sound  channels  and  for  the  wave  numbers  and  cutoff  frequencies 
associated  with  these  modes.  The  number  of  trapped  modes  is  shown  to  depend  on  the  gross 
features  of  the  sound  channel,  while  the  distribution  of  mode  wavenumbers  depends  to  a  greater 
degree  on  the  details  of  the  profile  shape.  Results  are  presented  for  the  square  (isovelocity- 
channel)  profile,  the  parabolic  profile,  and  the  bilinear  profile.  An  example  shows  how  the  simple 
formulae  can  be  applied  to  a  realistic  ocean  environment. 

While  the  analysis  is  presented  in  terms  of  underwater  acoustics,  the  results  are  applicable 
to  other  areas,  such  as  transmission  in  an  inhomogeneous  waveguide  or  to  the  solution  of  the 
Schrodinger  equation.  Most  of  the  results  presented  here  have  been  known  for  some  time2,  but 
what  is  naw  is  that  the  formulae  for  some  of  the  more  complicated  waveguides  can  be  put  in  the 
same  functional  form  as  the  well  known  formulae  for  the  ideal  waveguide.  The  physical 
parameters,  such  as  frequency,  depth,  and  sound  speed  are  easily  distinguished  from  the  details 
of  the  shape  of  the  sound  speed  profile,  which  can  be  treated  as  a  dimensionless  quantity.  The 
very  simple  form  of  the  expressions  makes  them  useful  for  back  of  the  envelope  calculations  or  for 
use  with  a  pocket  calculator.  Moreover,  using  the  same  functional  form  for  the  expressions  allows 
the  effect  of  the  shape  of  the  sound  speed  profile  to  be  easily  seen. 


2  The  normal  mode  equation 

The  normal  mode  equation  can  be  written  as, 

u„"(z)  *  Co)2/c2(z)-k2]un(z)  =  0  (1 ) 


where, 

z  is  the  depth  coordinate  (increasing  with  depth  from  the  surface), 

(4  a  27Tf  is  the  angular  frequency, 

c(z)  is  the  sound  speed  as  a  function  of  depth, 

k„  is  the  wave  number  or  eigenvalue, 

u„(z)  is  the  normal  mode  function, 

and  u"  denotes  the  second  derivative  of  u  with  respect  to  z. 
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In  general  one  wants  to  determine  the  normal  mode  wave  numbers  k„  end  the  associated 

mode  functions  u„(z),  subject  to  certain  boundary  conditions.  Notice  the  quantity  [U)2/c2(z)-k£] 
which  will  be  important  in  the  discussion  later;  in  particular,  it  is  equal  to  zero  at  a  turning  point, 
where  co/c(z)  «  it,,. 


SOUNO  SPEED 


Figure  1.  Normal  mode  solution  (heavy  dashed  line)  superimposed  on  a  sound  speed  profile 
(heavy  solid  line). 

Figure  1  shows  a  sound  speed  profile  (heavy  solid  line)  with  a  normal  mode  function  (heavy 
dashed  line)  superimposed  (with  arbitrary  amplitude)  at  the  appropriate  phase  velocity  vn  *  co/kn. 
Notice  a  number  of  things  about  the  mode  function: 

1)  At  the  turning  points  z,  and  z2,  0)/c(z,)  «0)/c(z2)  *  k„.  These  are  the 
classical  turning  points  for  tha  equivalent  ray  travelling  in  the  sound  channel. 

2)  At  the  air-water  interface  the  pressure,  and  hence  the  mode  function,  is 
zero;  i.e.  u„(0)  *  0. 

3)  Between  0  and  zt  the  normal  mode  function  has  an  increasing 
exponential  type  of  behaviour. 

4)  Between  z2  and  *  tha  solution  has  a  decreasing  exponential  type  of 
behaviour. 

5)  Between  z1  and  z2  there  are  three  zero  crossings  since  this  is  the 
fourth  mode;  (in  general  the  n-th  mode  will  have  n-1  zero  crossings). 

Furthermore,  note  tha  sinusoidal  behaviour  between  z,  and  z2;  note  also  that  u„  at  sin(7f/4)  at 
the  upper  turning  point  and  u„  a  sin(n-1/4)7r  at  the  lower  turning  point. 


3  The  WKB  method 


If  th«  sound  speed  profile  is  changing  slowly  with  respect  to  an  acoustic  wavelength,  the 
WKB  approximation1  allows  the  solution  to  be  written  in  terms  of  a  slowly  varying  amplitude  r(z) 
and  a  monotonicaily  increasing  phase  <£(2): 

un(z)  »  M  r(z)  sin[$(z)]  (2) 


where 


Kz)  *  [w2/e2(z)  -  k2]'1/4 


(3) 


0(z)  » 


[co2/c2(z‘)  -  k^1/2dz‘ 


z,  <  2  <  z2 


(4) 


and  M  is  a  normalization  constant.  The  integral  in  Eq.  (4)  will  be  referred  to  as  the  phase  integral. 
Note  several  points: 

1)  the  term  in  square  brackets  is  the  same  term  that  appeared  in  Eq.  (1). 

2)  near  a  turning  point  r(z)  is  singular. 

3)  0(z)  is  well  behaved,  however,  and  can  be  used  to  determine  the  WKB 
eigenvalues  k„. 

4)  0(z,)  3  5,  at  the  upper  turning  point. 

5)  <£(z2)  *  n7T-52  at  the  lower  turning  point. 

Figure  2  shows  how  the  phase  5,  at  the  upper  turning  point  depends  on  boundary  effects. 
If  the  surface  is  a  pressure  release  one,  the  pressure  is  zero,  un(0)  =  0  and  the  phase  =  0;  if 
the  surface  is  rigid,  the  normal  derivative  of  the  pressure  is  zero,  i.e.  un'(0)  3  0,  and  the  phase 
5,  *  TT/2.  At  a  turning  point  the  phase  5,  is  between  0  and  7T/2;  5,  3  tr/4  is  the  usual  choice1. 
The  same  comments  apply  to  the  phase  at  the  tower  turning  point. 

in  the  WKB  method  the  total  phase  change  between  the  turning  points  is  given  by: 


[go2/c2(z)  -  k51/2dz  +  5,  *  52  « 


rrff. 


(5) 
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PRESSURE  RELEASE 


RIGID 


TURNING  POINT 


8,*0  »,  *  %  0  <  8(  <  r/t 

Figure  2.  Effect  of  various  boundary  types  on  the  phase  5. 

where  51  and  S2  are  the  phases  at  the  upper  and  lower  turning  points.  With  the  usual  choice  of 
5,  *  52  *  7T/4  the  WKB  eigenvalue  equation  becomes 


[&>2/c2<2)  -  k£],/2dz 


•  (n-1/2)7r 


(6) 


0  c™  c’  C1Z)  0  Cm  C1Z)  °,  C|"  0(21 


(a)  (b)  (c) 


Figure  3.  Three  simple  sound  channel  shapes:  (a)  the  square  isovelocity-channel  profile, 
(b)  the  parabolic  sound  channel,  and  (c)  the  bilinear  sound  channel. 


One  wishes  to  solve  equation  (6)  for  k„  which  appears  explicitly  in  the  integrand  and 
implicitly  in  the  limits  z1  and  z2.  However,  for  certain  sound  speed  profiles  c(z)  the  integral  can  be 
evaluated  analytically  and  an  expression  obtained  for  kn.  Three  such  sound  speed  profiles  are 
shown  in  Figure  3  :  (a)  the  square  isovelocity-channel  profile,  (b)  the  parabolic  sound  channel,  and 

(c)  the  bilinear  sound  channel.  Note  that  it  is  actually  c'2  rather  than  c  which  is  parabolic  or  linear. 
Some  notation  is  also  introduced  at  this  point:  cm  =  the  minimum  sound  speed  in  the  channel,  c,  * 
the  maximum  sound  speed  in  the  channel,  and  h  *  z^-zu  »  the  maximum  vertical  extent  of  the  sound 
channel. 


4  The  number  of  trapped  modes 

The  WKB  eigenvalue  equation  can  be  solved  for  the  number  of  trapped  modes  N  if  the  phase 
integral  can  be  evaluated  either  analytically  or  numerically.  Using  n-»N,  z,  =  Zy,  z2  *  z/,  &  o>/c. 

and  rearranging  Eq.  (6)  gives  h  »  z^  -  z„  and  x  *  z/h: 


N  a 


1/2*(2hf/cm)[1-c2 


{[c2/c2(x)-c2/c2]/[1-c2/c2J>1/2  dx 


(7) 


equation  (7)  can  then  be  written  as 


N  ■  1  /2  ♦  (2hf/cw)  a  a 


(10) 


The  quantity  a  is  introduced  strictly  for  notationaf  convenience.  The  quantity  a,  however,  is 
related  to  the  shape  of  the  sound  speed  profile,  but  contains  none  of  the  physical  parameters 
such  as  the  frequency,  depth  or  sound  speeds.  Note  that  in  the  case  of  an  isovelocity  or  square 
profile  a  ■  1,  and  Eq.  (10)  gives  the  classical  formula  for  the  number  of  trapped  modes.  Figure  4 
graphically  illustrates  the  significance  of  the  integral  in  Eq.  (9),  where  the  vertical  extent  of  the 
channel  z^  -  z,  gets  mapped  into  the  range  0  to  1 ,  and  where  the  sound  speeds  between  cm  and  c. 
get  mapped  into  0  to  1  and  where  the  integrand  of  Eq.  (9)  (denoted  by  g(x)  )  is  enclosed  in  a 
square  box  of  unit  size.  The  integral  a  is  given  by  the  shaded  area,  which  can  be  calculated 
analytically  or  numerically  or  even  estimated  by  eye. 


? 
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C,  C(Z! 


g(x) 

WzM 


Z ‘  a* SHADED  AREA 

Figure  4.  The  sound  channel  and  the  associated  function  g(x);  the  shaded  area  a  is  defined 
in  Eq.  (9). 

5  Cutoff  frequencies 

The  cutoff  frequency  for  the  n-th  mode  can  be  obtained  by  rearranging  Eq.  (10)  to  give 


C*  *  [^(0-1/2)]  /[2haa] 


6  Wave  numbers 

For  the  three  sound  speed  profiles  shown  in  Figure  3,  the  eigenvalue  equation  (6)  can  be 
solved  for  the  wavenumbers  k„,  giving  results  of  the  form 


k2.u2/C2  -  A(n-1/2)P  (12) 

where  the  specific  values  of  A  and  p  are  given  in  Table  1.  For  the  three  profiles  considered  p 
varies  between  2/3  and  2,  while  the  corresponding  values  of  the  shape  parameter  a  varies  only 
between  2/3  and  1.  Moreover,  a  enters  equations  (10)  and  (11)  in  a  linear  fashion,  while  p 
appears  as  an  exponent  in  equation  (12).  Thus,  the  wave  numbers  are  much  more  sensitive  to  the 
profile  shape  than  are  the  trapped  modes  and  the  cutoff  frequencies. 
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Table  1 


Specific  values  of  a,  p  and  A  for  three  simple  profiles. 


Profile  Shape 

a 

P 

A 

Square 

1 

2 

(7T/h)2 

Parabolic 

1T/4  =  0.79 

1 

4acu/(cmh) 

Bilinear 

2/3 

2/3 

[37rw2a2/(2c2h)]2/3 

7  Example 

Table  2  shows  an  example  of  how  the  formulae  might  be  applied  to  a  realistic  ocean 
environment,  and  compares  the  wave  numbers  with  those  obtained  from  a  normal  mode  calculation. 
The  sound  speed  profile  approximates  a  typical  summer  sound  speed  profile  in  1 00  m  of  water  on 
the  Scotian  Shelf:  a  20  m  isovelocity  layer  of  speed  1 520  m/s  at  the  surface,  a  minimum  sound 
speed  of  1460  m/s  at  a  depth  of  40  m,  and  a  speed  of  1490  m/s  at  the  bottom.  The  table 
compares  the  wave  numbers,  or  phase  velocities,  obtained  using  the  bilinear  formula  with  those 
from  a  complete  normal  mode  calculation  at  200  Hz. 


Table  2 

Comparison  of  the  phase  velocities  obtained  from  equation  (12) 
with  those  from  a  complete  normal  mode  calculation. 


Mode  Number 


Phase  Velocities  (m/s) 

Normal  mode  Equation  (12)  Difference 


1  1465.72 

2  1471.34 

3  1475.99 

4  1480.13 

5  1484.94 

6  1487.48 


1465.46 

-0.26 

1471.43 

0.09 

1476.15 

0.16 

1480.29 

0.16 

1 484.09 

0.15 

1487.63 

0.15 

The  WK3  method  gives  a  good  approximation  to  the  normal  mode  wave  numbers;  four  or  five 
digits  accuracy  as  in  the  above  example  is  not  unreasonable.  In  fact  for  the  parabolic  profile  the 
WKB  and  exact  calculations  give  the  identical  results.  Provided  that  this  is  of  sufficient  accuracy, 
the  value  of  the  analytic  formula  is  obvious  from  a  computational  point  of  view. 


8  Discussion 

Equations  (10)-(12)  together  with  the  values  of  a,  p  and  A  given  in  Table  1  summarize  the 
results  of  this  paper:  simple  analytical  formulae  for  normal  mode  wave  numbers,  cutoff  frequencies, 
and  the  number  of  modes  trapped  in  a  sound  channel  of  simple  shape.  The  factor  of  1  / 2  appearing 
in  Eqs.  (10M12)  can  be  generalized  to  a  S  which  depends  on  the  boundaries  of  the  sound  channel 
as  well  as  the  type  of  turning  point  that  the  mode  "sees".  Equations  (10)  and  (1 1)  for  the  number 
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of  modes  and  the  cutoff  frequencies  are  useful  for  more  general  profiles  provided  that  the  integral 
of  Eq.  (9)  can  be  estimated. 

The  results  show  the  sensitivity  of  the  modes  to  the  shape  of  the  sound  speed  profile:  the 
number  of  modes  depends  on  the  shape  parameter  a  of  the  sound  speed  profile,  while  the 
distribution  of  mode  numbers  is  more  sensitive  to  the  details  of  the  profile  shape. 
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on  exact  and  WKB  solutions  are  presented  for  the  number  of  modes  trapped 
in  some  simple  sound  channels  and  for  the  wave  numbers  and  cutoff 
frequencies  associated  with  these  modes .  The  number  of  trapped  modes  is 
shown  to  depend  on  the  gross  features  of  the  sound  channel,  while  the 
distribution  of  modal  wave  numbers  depends  to  a  greater  degree  on  the 
details  of  the  sound  speed  profile  shape.  s 
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